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Abstract 

The method of QCD sum rules has proven to be particularly useful in heavy quark 
physics, where a small distance scale is provided by the inverse heavy quark mass. We 
present two new examples for the application of this method to B-physics: a calculation 
of the heavy-to-light form factors B — > ir, K and an estimate of the weak decay amplitude 
for B — > J/tpK beyond the factorization approximation. 
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1 Introduction 



Introduced fifteen years ago, the method of QCD sum rules |jj has become one of the most 
useful and reliable calculational schemes to derive hadronic properties from first principles 
with a minimal number of additional phenomenological assumptions. Originally, focusing on 
the operator product expansion (OPE) of two-point correlation functions one aimed at the 
determination of masses and coupling constants of ground states. The framework has then 
been extended and also applied to other dynamical characteristics of hadrons such as form 
factors and decay amplitudes. These extensions are based on OPE of three- and four-point 
correlation functions. Nevertheless, no new input parameters appear which is one of the most 
attractive features of this approach. All essential parameters such as vacuum condensates, 
coupling constants and thresholds of higher states are fixed by corresponding two-point sum 
rules. 

On the other hand, straightforward application of QCD sum rules in their original version 
to many-point correlation functions runs into a number of technical problems. First of all, 
the correlation functions are to be calculated in the deep Euclidean region. The result is then 
continued to physical momenta with the help of dispersion relations. In the presence of many 
channels, this continuation is not at all trivial. In some cases like for on mass-shell hadronic 
vertices and amplitudes one must rely on extrapolations which are plausible but not rigorous. 
Furthermore, the OPE for n-point correlation functions is usually represented by Feynman- 
type diagrams which are used to calculate the Wilson coefficients in perturbation theory. 
These diagrams grow rapidly in number with n and become rather complicated. Finally, in 
order to include higher resonances and continua of states into the sum rules, one often has to 
resort to models based on the concept of quark-hadron duality. These models usually work 
well for two-point correlation functions, but may be too crude already in the case of double 
dispersion integrals. 

For these and other reasons, it is very important to develop alternative, more economical 
methods which avoid these problems. Clearly, starting from vacuum correlation functions, 
applying OPE and deriving sum rules in terms of local vacuum condensates is only one way 
to account for nonperturbative quark-gluon dynamics. Other approaches leading to different 
versions of QCD sum rules include the external field method @-[|J and induced vacuum 
condensates, the concept of nonlocal condensates 0] , and the combination with light-cone wave 
functions 0-[[[(J. Each of these versions has advantages and disadvantages. The comparison 



of the results may shed some light on the accuracy and reliability of the different approaches. 

In this report we want to illustrate some of the current trends in the development of the 
QCD sum rule method by presenting two recent calculations relevant for B-physics. In section 
2 we discuss the calculation of the B — > 7r and B — > K form factors using an expansion in 
terms of pion and kaon wave functions on the light-cone with increasing twist, and compare 
the results with the predictions of conventional QCD sum rules. Then, in section 3, we study 
the problem how to calculate the amplitude of the prominent weak decay B — > J /ipK beyond 
the usual factorization approximation. In other words, we undertake the difficult task to 
estimate the nonfactorizable contributions to the decay amplitude originating from quark- 
gluon interactions at scales between the heavy b-quark mass and the hadronization scale. 
Following an idea suggested in ref. [11] for D meson nonleptonic decays, we try to extract 
this amplitude from a suitable four-point correlation function treated within the framework 
of conventional QCD sum rules. 
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2 QCD sum rules on the light-cone for the B — > 7r, X 
form factors 

The standard way to derive QCD sum rules for the transition amplitude < > between 

given ground state hadrons A and B starts from the vacuum correlation function of three 
currents, < 0|T{j^(x)j(0)j£(?/)}|0 >, where ]a (js) is chosen such as to create the hadron A 
(B) from the vacuum and T denotes the time-ordered product. A double dispersion relation in 
p\ and p 2 B , pa (pb) being the four-momentum in the A (S)-channel, is then used to relate this 
correlation function calculated in the Euclidean region with its imaginary part containing the 
desired transition amplitude. The remainder of the calculation is essentially the same as in the 
case of two-point QCD sum rules: derivation of OPE in terms of local quark-gluon condensates, 
subtraction of excited state and continuum contributions and Borel transformation. Since the 
quark-gluon condensates are universal, they can be taken from any QCD sum rule analysis. 
Furthermore, the coupling constants /a ~< 0|ja|^4 > and fs ~< 0\Jb\B > as well as the 
thresholds of higher states in the A and B channels are fixed by the two-point sum rules for 
the correlation functions < | ]a3a I > and < | ]b3b | >, respectively. Hence, as 
already pointed out, one has essentially no new parameters. This procedure has been applied 



to a variety of problems such as the pion form factor |12, 13], radiative charmonium transition 



rates |T4]]-[]l(|, and semileptonic form factors of heavy mesons []17l-|p2 



Here we demonstrate an alternative method which may be used in cases where one of the 
hadrons, A or B, is a light meson. An important example is provided by the B — ► 7r form 
factors f^r which determine the matrix element 

v TV 

< 7r\u lfl b\B >= 2/+(p 2 )g M + [f+(p 2 ) + /-(p 2 )K (1) 

with p + q, q and p being the B and 7r momenta and the momentum transfer, respectively. 
In the following, we concentrate on the phenomenologically more interesting form factor /+. 
The corresponding B — > K form factor can be treated in parallel by obvious formal 
replacements. Numerical results will be shown for both /+ and f^. 

Instead of investigating the vacuum averaged correlation of the b — > u transition current 
with two other currents carrying the B and ir quantum numbers, we consider the matrix 
element 

F„(p, q)=i [ d A xe ipx < 7r(q) | T{u(x) 7lt b(x),b(0)vy B d(0)} | > (2) 



= F{{p + q)\p 2 )q ll + F{{p + q)\p 2 )p ll 

between the vacuum and an on-shell pion state. This object is represented diagrammatically 
in Fig. 1. The pion momentum squared, q 2 = m 2 , vanishes in the chiral limit adopted 
throughout this discussion. Moreover, the light u and d quarks forming the pion eventually 
propagate to large distances. In contrast, the b-quark still propagates far off-shell provided 
that (p + q) 2 is taken sufficiently large and negative, and the time-like momentum transfer 
squared p 2 is far from the kinematical limit, p 2 = m 2 B . 

Formally, contracting the 6-quark fields in (^) and keeping only the lowest order term, i.e. 
the free 6-quark propagator, yields 

Fto + M-if J ge*"^ ? , (3, 
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where 

<P a {x\ q-x)=< ir(q) \ u(x)T a d(0) | > , (4) 

T a denoting certain combinations of Dirac matrices. This approximation corresponds to Fig. 
la. The high virtuality of the 6-quark propagating between the points x and guarantees 
that the u and d quarks are emitted at almost light-like distances. In that case, it is justified 
to keep only the first few terms in the expansion of the matrix elements (f|) around x 2 = 0, 
that is near the light-cone: 



[x 



2 q-x) = J2€(q-x)(x 2 r . (5) 



Logarithms in x 2 which may also appear in (||) are disregarded for simplicity. These terms can 
be consistently treated by means of QCD perturbation theory. They give rise to normalization 
scale dependence. Because of translation invariance, the coefficients 0™ must have the form 

0a(<7 ' x ) ~ / du(pa(u)exp(iuq ■ x) . (6) 
Jo 

Inserting (H) and (^) into (^) and integrating over x and P, one obtains, schematically, 

It is thus possible to calculate the invariant function F with reasonable accuracy in the kine- 
matical region of highly virtual b-quarks provided one knows the distribution functions 
at least for low values of n. These distributions contain essential information about the dy- 
namics at large distances and play a similar role as the vacuum condensates in the conventional 
OPE of vacuum correlation functions. 

As it turns out, the distributions are nothing but light-cone wave functions of the 

pion the twist of which is growing with n. These wave functions were introduced long ago in 



the context of hard exclusive processes [p3||-||26|. The leading twist 2 wave function is defined 
by 

< 7r(g) | u{x)j^ 5 Pexp{i f da x lM A tl {ax)}d{^) | >= -if^ f due iuq ' x y^{u) , (8) 

io J 

where the exponential factor involving the gluon field is necessary for gauge invariance. The 
asymptotic form of 4> n (u) can be obtained from perturbative considerations and is well known 

<Pn(u) = 6u(l — u). (9) 



Over the years a great deal has been learned about these wave functions |2"?|-|3"2]]. All im- 
portant twist two, three and four wave functions have been identified and their asymptotic 
form has been determined. Besides two-particle wave functions reflecting the quark- ant iquark 
structure of mesons, also three-particle wave functions associated with the quark- ant iquark- 
gluon component have been studied. Very important for practical applications, nonasymp- 
totic corrections have been estimated considering expansions in orthogonal polynomials and 
renormalizing the coefficients of this expansions in QCD perturbation theory. Furthermore, 



conventional two-point QCD sum rules have been used P5 to fix the first few moments of 



the low twist wave functions. Based on these estimates and improvements, various models 
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for the nonasymptotic form of these functions have been suggested ( see, for example, refs. 
[10,28,32]). 

The approach outlined above together with the accumulated knowledge of pion and kaon 
wave functions has been employed in ref. [6] to calculate the B — > tt form factor at zero 
momentum transfer, in ref. [7] to get the D — > tt form factor, and in ref. [8] to derive the 
B — > 7i and B — > K form factors including the momentum dependence. 

In our calculation of /+ and /j£ we have included quark-antiquark wave functions up to 
twist four. In addition, we have also estimated the first-order correction to the free 6-quark 
propagation shown in Fig. lb which involves quark- ant iquark-gluon wave functions. On the 
other hand, the perturbative 0(a s ) corrections indicated in Figs, lc and Id have not been 
evaluated directly, but have been taken into account indirectly as explained below. The direct 
calculation of these corrections is on its way. 

In order to extract the desired form factor from the resulting invariant function F((p + 
q) 2 ,p 2 ) sketched in (^) we employ a QCD sum rule with respect to the I?-meson channel. 
Writing a dispersion relation in (p + q) 2 , we approximate the hadronic spectral function in 
the B-channel by the pole contribution of the B meson and a continuum contribution. In 
accordance with quark-hadron duality, the latter is identified with the spectral function derived 
from the QCD representation (0) above the threshold (p + q) 2 = Sq. Formally, subtraction of 
the continuum then amounts to simply changing the lower integration boundary in (7) from 
to A = (m 2 — p 2 )/(so — p 2 ). After Borel transformation one arrives at a sum rule for the 
product f^fs where f B is the B meson decay constant. It is important to note that the 
numerical values to be substituted for m b , f B and the threshold s are interrelated by the 
QCD sum rule for the correlation function < | T {b{x)^^u{x) , ^(0)756(0)} | >. We stress 
that the latter sum rule should be used without 0(a s ) corrections in order to be consistent 
with the neglect of these corrections in the sum rule for f^f B (see also ref. [20]). 

The final expression of the form factor /+ is given by 



2f B m 2 B JA u M 2 uM 2 1 

Mm) + |.„ Mu) + ^(„)( 2+ !^)], (!0) 

where M 2 is the Borel parameter and /1 = m^./(m u + md). While in the concrete calculation of 
ref. [8] the twist 2 wave function ip n is corrected for nonasymptotic effects p8|,p2|, the twist 



3 wave functions (p p and (p ff are taken in their asymptotic form. Numerically, the higher twist 
contributions turn out to be more important than the nonasymptotic corrections to the leading 
twist wave function. For brevity, in fllO|) we have omitted the contributions from the twist four 
quark-antiquark wave functions and the term involving the twist three quark-antiquark-gluon 
wave function. These contributions are quantitatively negligible. Further details can be found 
in ref. [8]. 

Numerical results for f^(p 2 ) are plotted in Fig. 2a for two values of the Borel parameter 
M characterizing the fiducial range of the method. As can be seen, the predictions are quite 
stable under variation of M as long as p 2 is not getting too close to m 2 B . For comparison, 
we also show results obtained by other methods. Furthermore, it is interesting to note that 
fn{p 2 ) is not very sensitive to the precise shape of the leading twist wave function (pn(u) at 
least at p 2 < lOGeV 2 . We have checked this by replacing the nonasymptotic two-humped 
wave function |23| used for Fig. 2a by the simple asymptotic wave function given in (PI). 
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The B — > K form factor /j£ was calculated analogous to /+ using the leading twist wave 
function 99 ^ from ref. [28]. The higher twist wave functions are left unchanged. The numerical 
result is presented in Fig. 2b. Although is not accessible directly in semileptonic B decays, 
it determines the factorizable part of the nonleptonic decay amplitude for B — > J/ipK and 
is therefore phenomenologically very important. Later, in section 3 we shall make use of the 
value 

/+( m 2) = 0.55 ±0.05. (11) 

In conclusion, we emphasize that light-cone sum rules such as the ones exemplified in this 
section represent a well defined alternative to the conventional QCD sum rule method. In this 
variant, the nonperturbative aspects are described by a set of wave functions on the light-cone 
with varying twist and quark-gluon multiplicity. These universal functions can be studied in a 
variety of processes involving the it and K meson, or other light mesons. The most important 
advantage of the light-cone sum rules is the possibility to take hadrons on mass-shell from the 
very beginning. One can thus avoid the notorious model-dependence of extrapolations from 
Euclidean to physical momenta in light channels. Furthermore, in many cases the light-cone 
approach is technically much easier than the conventional QCD sum rule technique. Finally, 
the light-cone method is rather versatile. It can also profitably be employed to calculate 
heavy-to-light form factors such as B — ► p and B — ► K*, amplitudes of rare decays such 
as B —>■ A'*7 0, and hadronic couplings such as B*Bn [Q. Just in passing we mention 
that the B*Btt coupling can be extracted from the correlation function (2) by performing a 
second Borel transformation in p 2 . Needless to say, light-cone sum rule are equally useful in 
calculating properties of D mesons and also of light hadrons [jT0|l , [j32"l . 

The main problem to be solved if one wants to fully exploit the light-cone approach is 
a reliable determination of the nonasymptotic effects in the wave functions. In this respect, 
measurements of hadronic form factors, couplings etc. can provide important information. 
A second, mainly technical problem, concerns higher order perturbative corrections, such as 
those shown in Figs, lc and Id, and the possible occurrence of large logarithms. 



Weak decay amplitude for B —> J/i/jK beyond fac- 
torization 



Since the earliest estimates El [BH], nonleptonic two-body decays of heavy mesons are usually 



calculated by splitting the appropriate matrix element of the weak Hamiltonian into a product 
of a semileptonic transition form factor and a decay constant of one of the final mesons. On the 
theoretical side, one still lacks a strict proof of this recipe. Empirically, with the accumulating 



data on D decays, it soon became clear that naive factorization fails ||37|| . In order to achieve 
agreement with experiment it is necessary to let the coefficients of the OPE of the weak 
Hamiltonian deviate from their values predicted in short-distance QCD. Phenomenologically 



|38| 1, the coefficients are reinterpreted and treated as free parameters to be determined from 



experiment. For D decays, it has been shown |59j that the fitted values can be reconciled 
with the short- distance expectations if the matrix elements of the local operators appearing 
in the OPE are expanded in 1/N C and if only the leading terms are kept. This observation 
gave rise to the rule of discarding nonleading in 1/N C terms. Factorization combined with this 
rule leads to a satisfactory description of the majority of two-body .D-decays. 

For B decays, an instructive example is provided by the illustrious decay mode B — > J/ipK. 
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The effective weak Hamiltonian responsible for this decay is given by 

H w = -^^^{dOi + c 2 2 } , (12) 

O x = (^c)(cl>), 2 = (cTPc)(sT p b) (13) 

where G is the Fermi constant , V cb and V cs are the relevant CKM matrix elements and 
Tp = 7p(l -75)- The coefficients ci(/i) and c 2 (/z) incorporate the short-distance effects arising 
from the renormalization of Hw from fi = my/ to fi = 0(m b ). These coefficients are known 
to next-to-leading order (NLO) [[4(]] , j41| . While 2 already possesses the appropriate flavour 
structure, this is the case for 0\ only after Fierz transformation: 

1 = l -0 2 + 20 2 , (14) 

where 

d 2 = (cT^c)(sT p ^b)} . (15) 
Assuming factorization, the matrix element of H\y for B — > J/ipK simplifies to 

< Jji)K I (c 2 + ^)0 2 + 2 Cl 2 I B > 

= (c 2 + — ) < J/V> I cT p c I >< K I sY p b I B > . (16) 
3 

Note that in this approximation 

< J/ipK \6 2 \B >=< J/ip I cr^yc I >< K I sY p ^-b I B >= (17) 

because of color conservation. Moreover, using a parametrization of < K \ s'jpb \ B > 
analogous to (||) and < J/ijj \ cy p c | >= m^f^ep one readily obtains the following expression 
for the decay amplitude: 

A(B - J/^K) = V2GV cb V; s (c 2 + j)Uf£(p 2 = m%)m^e ■ q) . (18) 

The corresponding branching ratio is given by 



BR{B -> J/^K) = —\V cb V;f(c 2 + ^fflft\v 2 = 



x m|(l - ) 3 r B (19) 

where r B is the B meson lifetime. 

Taking for the B — > K form factor the value (|TT|) from ref. [8] which agrees with other 
calculations, and using = 409 MeV, V cb = 0.04, and = (1.489 ± 0.038) ps we obtain 

0.005% < BR(B -> J/^fsT) < 0.013% . (20) 

The range of this prediction mainly reflects the theoretical uncertainties in the coefficients C\ 
and c 2 at /i = m b , given in ref. [42]: 

{d = 1.115, c 2 = -0.255} {d = 1.146, c 2 = -0.312} . (21) 



6 



As a matter of fact, the combination ci + ci/3 turns out to be particularly sensitive to the 
precise values of C\ and c 2 , as pointed out in ref. [43]. Despite of this uncertainty, the 
branching ratio expected from naive factorization is undoubtedly much lower than the recent 
experimental results [fHf : 



BR(B~ -> J/ij}K-) = (0.11 ± 0.015 ± 0.009)% 

BR(B° -> J/^AT°) = (0.075 ± 0.024 ± 0.008)% . (22) 

On the other hand, applying the rule of discarding the nonleading in 1/N C terms, that is 
dropping the term ci/3 in (0), yields 

0.065% < BR(B -> J/V>#) < 0.095% (23) 

in good agreement with (|22|). 

Yet, from this observation one cannot conclude that the l/iV c -rule generally works in B 
decays. The CLEO analysis [fHfl , for example, shows that factorization is consistent with the 



measurements provided the short distance coefficients c\ + C2/3 and C2 + Ci/3 are replaced 
by two a priori unknown parameters a± and 02, respectively. Similarly as in D decays, these 
coefficients are found to be universal, at least on the basis of the two-body B decays analysed 
so far. However, unlike in D decays the data favour a value 02 > which is inconsistent with 
the l/iV c -rule implying a 2 — c 2 < 0. ( see also the discussion in ref. [45]). 

The present situation is very unsatisfactory and rather confusing. Theoretically, there 
seems to be no trustworthy reason why factorization should work. Moreover, on dynamical 
grounds j|8] the universality of a% and 02 is rather surprising than expected. Also, the empirical 



facts do certainly not prove factorization. Just the contrary may be the case. For example, in 



the matrix element (|16 ) the factorizable term proportional to ci/3 may cancel with sizeable 
nonfactorizable contributions neglected in eqs.(^) to (^0|) . Effectively, this would lead to the 
result (p3|) and explain the success of the 1/N C rule. Such a cancellation was first advocated in 
ref. [38] and then shown in ref. [11] to actually take place in two-body D-decays by estimating 
the nonfactorizable contributions to the decay amplitudes using QCD sum rule techniques ( 
for a more recent discussion see refs. [46,47]). This cancellation may be different in B decays. 
It may occur in certain channels and not in others. The nonfactorizable contribution may also 
happen to overcompensate the factorizable amplitudes nonleading in 1/N C , and thus change 
even the sign of the effective coefficient a<i- Clearly, one has to explain the coefficients a\ and 
02 before one can claim some theoretical understanding. 

Recently, we have started an attempt in this direction f4"8fl . We are investigating the 
problem of factorization in B decays using B —>■ J / ipK as a study case. Here we briefly explain 
our approach and present some preliminary results. Following the general idea suggested in 
ref. [11], we consider the four-point correlation function : 

lWff,p) = / d'xd'yd'ze^y < I T{j«(x)jt(y)H w (z)j£(0)} | > , (24) 

where j% — u^y^s , jt = crf v c and jf = bi^ 5 u are the generating currents of the mesons 
involved and Hy/ is the effective weak Hamiltonian given in (|i"2|). The four-momenta assigned 
to the K , J I ip and B channel are q, p and P = q + p, respectively. The momentum transfer 
in the weak interaction point z is zero. The bare diagram associated with (p4[) is depicted in 
Fig. 3a. 
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The correlation function is then calculated by means of the QCD operator product 
expansion. The appropriate Euclidean region located far enough from the physical thresholds 
in all three channels is q 2 < —lGeV 2 , p 2 < and P 2 < 0. For the J/ip and B channels, 
the conditions are more relaxed than for the .fT-channel, since the former are protected by the 
large c and b quark masses. The light quarks are considered massless. In the OPE for ([24] ) 
we include all local operators up to dimension six. The Wilson coefficients are calculated to 
lowest nonvanishing order in a s . 

In this approximation the nonfactorizable contributions to (f24[) only arise from the operator 
O2 given in (|T5|). This is just the operator the matrix element of which vanishes by factorization 



as pointed out in (|17| ). However, nonperturbative gluon exchange between the loop and the 
triangle indicated in Fig. 4 breaks factorization and gives rise to a finite contribution from O2. 
Fig. 4a shows one of the diagrams associated with the d = 4 gluon operator G° llv G ailv . Typical 
diagrams for the d = 5 quark-gluon operator qgG® LU \ a <j^ u q and the d = 6 four-quark operator 
are depicted in Figs. 4b and 4c, respectively. Additional d = 6 contributions corresponding 
to the diagram Fig. 4b emerge from the operators qV a qG^ x and qD a G^ x q, where V a and V a 
denote the proper covariant derivatives. 

Of course, there are also perturbative gluon corrections to the bare diagram, Fig. 3a, 
associated with the unity operator in the OPE. These are exemplified in Fig. 3b and 3c. 
The analogous corrections to two-point functions play an essential role in the corresponding 
sum rules. Therefore, they should also be included in a complete treatment of the four-point 
function (|24]) . Some of the perturbative gluon corrections (Fig. 3b) are in fact included, 
summed up in the coefficients ci(mfe) and C2(mt,) of Hw, namely those originating at short- 
distances between I /raw and 1/mb- However, there are other contributions (Fig. 3c) coming 
from larger distances, say, between 1/mb and l/m c or even 0(lGe^ _1 ), which are definitely 
not included in Hw- These may give rise to logarithms of the kind ln(mb/m c ) and induce 
important nonfactorizable contributions to the decay amplitudes. A systematic inclusion of 
all perturbative gluon effects requires three-loop calculations with massive quarks, a task 
which is postponed to later developments. Restricting in the first step of the analysis to the 
nonperturbative interactions sketched in Fig. 4, the correlation function (|24]) receives the 



following nonfactorizable contributions: 

n?? D = n£, G + + nJ^ G + + n*« , (25) 



where the tilde is a reminder to replace Hw in ([24]) by O2. The individual terms in ( p5|) have 
been derived from the corresponding diagrams of Fig. 4 in the form of Feynman integrals. 
Their explicit expressions can be found in ref. [48]. 

With the result (|25|) at hand, one can now proceed in constructing sum rules which will 
allow to determine the matrix element ( |TTD entering the correlation function through the 
resonance contribution 



ttb^wk = . < I 3% I K >< I g I JAA >< KJjj, I Q 2 I BxB \ j B | > 
^ ' {m 2 K - q 2 ){m% - p 2 )(m| - P 2 ) 



However, this problem is much more difficult than it sounds. One reason is that apart from 
the ground state contribution (|26| ) and analogous contributions from excited states, one also 
has to take into account contributions to the spectral functions from continuum states with 
flavour quantum numbers which differ from those of the respective currents. Such intermediate 
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states emerge in every channel of the correlation function due to final state interaction. For 
example, in the 5-meson channel the following contribution appears inevitably: 



< | jfc | K >< | jt | J/ip >< KJ/if) | n DD* s n >< v DD* s n \ 6 2 j B | > 



" «-g 2 )K-p2 )(m ^_ p2) I ) 

Here, the questionable intermediate state carries the quantum numbers of a virtual DD* 
state. It is created by weak interaction and converted into the J/ipK final state by strong 
interaction. The problem is that this continuum of states has a mass-threshold below the 
5-meson pole. It can therefore not be suppressed by Borel transformation and subtracted 
away similarly as outlined in section 2 for the normal continuum in the 5-channel starting 
at a threshold So > m%. On the contrary, this unwanted contributions will be enhanced 
by Borel transformation relative to the contribution from the ground state B roughly by a 
factor exp[(m 2 — Am 2 c )/M 2 ] which is quite large at characteristic values of the Borel parameter 
M 2 ~ [m\ — m 2 h ). 

Such dangerous contributions are also present in D decays. They have been estimated in 
ref. [11] using a simple model and found to be unimportant. Unfortunately, in the present 
case the solution of this problem seems to be less trivial. We have identified the analogous 
contributions in the QCD diagrams of Fig. 4 and find that they are numerically very important 
after Borel transformation. Indeed, examination of the analytical properties of the diagrams 
in Fig. 4 shows that in addition to a pole at P 2 = ml the four quark condensate diagram, Fig. 
4c, has a discontinuity at P 2 > Am 2 c connected with the four-quark intermediate state with the 
flavour combination usee. This is just the flavour composition of the virtual DD* continuum 
in Q2"7p. It seems to be reasonable to again use the quark-hadron duality principle in order 
to cancel this piece of (|25|) against the unwanted hadronic contribution (p7|) . Practically, this 
results in neglecting the four-quark condensate diagram altogether. 

From here on, we follow the usual procedure. We perform a Borel transformation in P 2 
and take moments in p 2 as is usually done in the charmonium channel. The momentum 
squared in the i^-meson channel is kept fixed at Q 2 = —q 2 ~ lGeV 2 . In order to take into 
account the higher states in the resulting sum rules we adopt for simplicity a two-resonance 
description in each of three channels, i.e. above the effective thresholds soi we describe the 
hadronic spectral functions by one effective resonance with a mass equal to a/s^. Normalizing 
the matrix element ( |T7| ) so that the amplitude ( |T^| ) is proportional to combination 

a 2 f+ = (c 2 + ^)f+ + 2cj, (28) 

we obtain for the nonfactorizable contribution 

/= -(0.045 -0.075) . (29) 

The range reflects our estimate of the uncertainty in the QCD calculation and the continuum 
subtraction. 

Our analysis shows that the factorizable, nonleading in 1/N C term and the nonfactoriz- 
able term in (|28|) are opposite in sign. Quantitatively, the nonfactorizable matrix element 
is considerably smaller than the factorizable one, — 0(10%). Nevertheless, the 

nonfactorizable contribution to the decay amplitude is very important because of the large 
coefficient, 2c\j{c 2 + ci/3) ~ 20-j-30. In fact, if |/| is close to the upper end of the range given 
in (|29|), the nonfactorizable contribution almost cancels the factorizable one proportional to 
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Ci/3, thus leading to the prediction ( p3]) which agrees with experiment. This is exactly the 
scenario anticipated by employing the l/N c -rvle. 

At the same time, in our approach the nonfactorizable contributions are expected to be 
channel-dependent on quite general grounds. This expectation is corroborated by our prelimi- 
nary analysis of other nonleptonic 5-decays. In other words, as is obvious from eq. fl28|) there 
is no theoretical reason to expect a single universal value for the effective coefficient a%. Even 
the sign of 02, depending on the sign of / in a given channel, may vary. Universality can at 
most be anticipated for certain classes of decay modes, such as B — > Dn or B — > DD, etc. 
Moreover, there is no simple relation between B and D decays in our approach since the OPE 
for the corresponding correlation functions significantly differ in the relevant diagrams and in 
the hierarchy of mass scales. We hope to be able to clarify this issue further. Theory seems to 
predict a much richer pattern in two-body weak decays than what is revealed by the present 
phenomenological approach to the data. 

4 Conclusion 

We have discussed applications of QCD sum rules in different variants to form factors and 
exclusive decay amplitudes of B mesons, both being important issues in I?-physics. There 
are many other interesting aspects which can be studied in the framework of QCD sum rules 
and which have not been touched upon here. We just mention the dependence of form factors 
and couplings on the heavy quark mass. Here, the sum rule technique can be employed 
complementarily to the heavy-quark effective theory. From a theoretical point of view, heavy- 
light bound states such as the B meson are very suitable systems for studying the dynamics 
of the light quark and gluon degrees of freedom and for probing nonperturbative methods 
in QCD. It will be interesting to see whether or not the QCD sum rule technique passes 
these tests as successfully as the various tests in the past. For the time being, it is the only 
viable approach to some problems such as the calculation of weak decay amplitudes beyond 
factorization. Even the capabilities of current lattice calculations do not allow to solve this 
problem (as remarked in ref. [49]). 
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Figure Captions 



Fig. 1: QCD diagrams contributing to the matrix element (2) involving (a) quark- ant iquark 
light-cone wave functions; (b) three-particle quark-antiquark-gluon wave functions; (c) and 
(d) perturbative 0(a s ) corrections. Solid lines represent quarks, dashed lines gluons, wavy 
lines are external currents. 

Fig. 2: _B-meson form factors calculated from light-cone sum rules: (a) the form factor 
fn(p 2 ) of B — > 7r transitions and (b) the form factor f^ip 2 ) of B — > K transitions at M 2 = 
lOGeV 2 (upper solid curves) and M 2 = lbGeV 2 (lower solid curves). The quark model 
predictions from ref. [38] (dash-dotted curves) and the QCD sum rule results for f^{p 2 ) from 
ref. [20] (dashed curve) and for f^(0) from ref. [6] (arrow) are shown for comparison. 

Fig. 3: Diagrams associated with the correlation function (24): (a) bare diagram with 
H w replaced by 2 ; (b) and (c) diagrams corresponding to perturbative gluon corrections. 

Fig. 4: Diagrams associated with (a) the gluon condensate, (b) the quark-gluon conden- 
sate and (c) the four-quark condensate contributions to the correlation function (24) with Hw 
replaced by 2 . 
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